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ABSTRACT 

Context. MWC 656 is the recently discovered first binary system case composed of a Be-type star and an accreting black hole. Its 
low X-ray luminosity indicates that the system is in a quiescent X-ray state. 

Aims. The aim of our investigation is to establish if the MWC 656 system has detectable radio emission and if the radio characteristics 
are consistent with those of quiescent black hole systems. 

Methods. We used three archived VLA data sets, one hour each at 3 GHz, and seven new VLA observations, two hours each at 
10 GHz, to produce very high sensitivity images down to ~1 pjy. 

Results. We detected the source twice in the new observations: in the first VLA run, at apastron passage, with a flux density of 
14.2 ± 2.9 yuJy and by combining all together the other six VLA runs, with a flux density of 3.7 ±1.4 /jJy. The resulting combined map 
of the archived observations has the sensitivity of Icr = 6.6pJy, but no radio emission is detected there. 

Conclusions. The radio and X-ray luminosities agree with the behaviour of accreting binary black holes in the hard and quiescent 
state. In particular, MWC 656 in the Lx, Lr plane occupies the same region as A0620-00 and XTE J1118±480, the faintest known 
black holes up to now. 
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1. Introduction 

In 2010 the gamma-ray source AGL J2241+4454 was de¬ 
tected by the AGILE satellite ( |Lucarelli et al.||2010[). Within 
the positional error circle of the satellite Williams et al. ( 2010| l 
found the emission-line Be star MWC 656 (also known as HD 
215227). The first suggested classification for MWC 656 was 
B3 IVne+sh, where n and e indicate broad lines and Balmer 
emission, and sh denotes the presence of shell lines, indicated 
by a sharp central absorption in the line. The V/R ratio of 
the Hy emission changed significantly in one day, which is 
unusually fast for most Be stars and Williams et al. ( 2010| l 
pointed out some resemblance of the spectral properties of this 
newly identified gamma-ray source to the gamma-ray binary 
LS I +61°303 and its Ha variations. The star MWC 656 is at a 
Galactic latitude ofb- -12°; at a distance of d - 2.6 + 0.6 kpc 
the star results well below the Galactic plane at z = -560 + 200 
pc. This is an extreme distance for a normal OB star and 
I Williams et al.| ( |2010| l suggested a runaway star, which reached 
the current position by a supernova explosion of the compan¬ 
ion star. Timing analysis of variations in the optical photom¬ 
etry, interpreted as orbital modulation, resulted in a period of 
60.37±0.04 days with the epoch of maximum brightness at HJD 
2453243.3± 1.8. 

The binary hypothesis was tested with a radial velocity 
study by Casares et al.| (|2012|l. R adial velocities folded with 
the 60.37 d of Williams et al.| ( |2010| ) revealed a sine-like 
modulation. Casares et al. (20121 by the rotational broadening 
V sin i ~ 346 km s“^ estimated the inclination of the orbit i >66° 
and a mass of the compact object in the range ~ 2.7 - 5.5 
Mq. Moreover, Casares et al. (2012 1 found that the main param¬ 
eters of the Ha emission line (equivalent width, full width at half 


maximum, and centroid velocity) result to be modulated by the 
proposed orbital period. 

The turning point in this research came with the new opti¬ 
cal observations by Casares et al.| (2014i. To improve the radial 
velocity curve of the Be star, Eell emission lines from the in¬ 
nermost region of the Be disk were used instead of stellar broad 
absorption lines contaminated by the circumstellar wind emis¬ 
sion lines. On the other hand, the spectra show an Hell 4686 A 
double-peaked emission line, indicating a disk, whose centroid 
is modulated with the same 60.37 day period, but in antiphase 
with respect to the radial velocity of the Be star, as obtained from 
the Eell line (their Eig. 3). This important observation therefore 
hints at an accretion disk around an invisible companion. The 
orbital solution resulted in Operiastron =0.01+0.10 (phase zero set 
to HJD 2453243.7). With new observations, the authors were 
also able to obtain a better spectral classification for the Be star, 
and determined a B1.5-B2 III classification. Given the mass of 
this star, 10-16 Mq , the implied companion mass is 3.8-6.9 Mq. 
This makes MWC 656 the first clear case of a Be-type star with 
a black hole companion. The case of LS I +61°303 is, in fact, 
still unclear because of the uncertainties in the values of mass 
function and inclination of the orbit (|Massi|2004l |Casares et al.| 
|2005l ). 

Observations with XMM-Newton at O = 0.08 by |Munar-| 
Adrover et al. ( 2014[ l proved that MWC 656 is indeed an X-ray 
binary system. The X-ray luminosity of ~ 10^'erg s ' points to a 
stellar mass black hole in quiescence. The quiescent X-ray state 
is one of the X-ray states of an accreting black hole similar to the 
low-hard X-ray state. Eor a compact object of a few solar masses, 
the low/hard state corresponds to Ly ~ 10^^ erg s ' but rnay dr op 
to Lx = 10^" ^ - 10^^ ^ erg s ' |McClintock & Remillard 20061 at 
its lowest phase, called the quiescent state. The X-ray luminosity 
in the quiescent and low/hard X-ray states correspond to a radia- 
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Table 1. MWC 656 observations and final parameters of maps. 


Epoch 

Date 

Julian Date 

O 

Synthesized beam 

rms noise 

Peak 

Flux±fr'' 

(day.month hh:mm) 



(0maj["] X P.A.[°]) 

(pjy beam^') 

(juJy bm“') 

(ply) 

X-band observations in 2015 

1 

22.02 19:15 

2457076.302 

0.49 

1.17 X 1.09; 101.7 

2.9 

9.9 

14.2 ±2.9 

2 

28.02 19:30 

2457082.313 

0.58 

1.15 X 1.08; 104.3 

3.8 


<11.4 

3 

07.03 22:10 

2457089.424 

0.70 

2.69 X 0.36; 58.8 

2.3 


<6.9 

4 

08.03 23:04 

2457090.461 

0.72 

1.62x 1.06; 84.6 

2.4 


<7.2 

5 

03.04 19:36 

2457116.317 

0.15 

1.89 X 1.09; 89.3 

2.2 


<6.6 

6 

10.04 17:01 

2457123.209 

0.26 

1.30x0.97; -20.6 

2.7 


<8.1 

7 

12.04 20:25 

2457125.351 

0.30 

1.10 X 1.02; -21.9 

2.3 


<6.9 

2-7 




0.84x0.70; -75.1 

1.02 

3.47 

3.7 ± 1.4 

1-7 




0.82x0.70; -75.1 

0.95 

4.15 

4.5 ± 1.2 

S-band observations in 2012 

1 

05.10 04:38 

2456205.693 

0.06 

0.48x0.41; 12.3 

12.6 


<37.8 

2 

15.10 08:03 

2456215.835 

0.23 

0.66 X 0.44; 89.6 

10.1 


<30.3 

3 

06.12 04:54 

2456267.704 

0.09 

0.61 X 0.50; 83.0 

11.1 


<33.3 

1-3 




0.51 X 0.45; 80.5 

6.6 


<19.8 

aer limits at three times the noise level. 






inefficient. 

“jet-dominated” accretion mode (Fender et al. pier and 32 different spectral windows (with bandwidths 


2003). In this mode, only a negligible fraction of the binding 
energy of the accreting gas is directly converted into radiation. 
Most of the accretion power emerges in kinetic form, as shown 


for Cygnus X-1 by Gallo et al. (2005 1 , and for AGNs through 


the relationship between the Bondi power and the kinetic lumi¬ 
nosity ( [Merloni & Heinz|2007[ and references therein). That is, 
during the low-hard and quiescent states the liberated energy of 
the accretion is thought to be converted into magnetic energy, 
which powers the relativistic jet observed in these states (jGalloj 
jet al.|2003[|Fender et al.|2004l|Gallo et al.|2006||Smith|2012|l. 

The origin of the X-ray emission during the low-hard and 
quiescent states is still controversial. The emission may be due to 
a Comptonizing corona and/or to the jet (see [Gallo et al.|2006| l. 
Nevertheless, it is well established that the X-ray emission is 
related to the radio emission from the jet (Gallo et al. 2003 1 2006 


[Merloni et al.|200^ . This important nonlinear scaling between 
X-ray and radio luminosities. Lx, L^, has been demonstrated to 
hold, with the addition of a mass term, across the entire black 
hole mass spectrum, from microquasars to AGN ([Merloni et al.[ 
[2D^ . 

Therefore, if a jet is associated with MWC 656, its radio 
emission should be consistent with the Lx, Lr relationship. The 
aim of our investigation is to establish if the system MWC 656 
has an associated radio emission and if the radio emission ful¬ 
fills the Lx, Lr relationship. In this letter, we present new radio 
observations in the direction of this system. Section 2 describes 
the observations and the data reduction. In Sect. 3 we report on 
our results. Section 4 describes MWC 656 in the context of the 
Lx, Lr relationship, and finally. Sect. 5 gives our conclusions. 


2. Observations 

We obtained seven new X-band (8 to 12 GHz) observations with 
the Karl G. Jansky Very Large Array (VLA) of the National 
Radio Astronomy Observatory (NRAO) in its B configuration. 
These observations were made under project 15A-013. The re¬ 
ceiver was used in semicontinuum mode with the 3 bit sam- 


MHz) were recorded simultaneously to cover the full band. 

Each observing session ran for two hours and was organized 
as follows. First, we spent 4.5 minutes on scans for instrument 
setups as recommended by NRACj^ Then, we observed a 5.5 
minute scan on the phase calibrator J2255-1-4202, and the large 
scan was performed to take the slewing time into account. Then 
we observed nine cycles of 10.6 minutes on the target followed 
by 1.0 minutes on the phase calibrator. We finished the observa¬ 
tion with a 5.5 minutes scan in the flux calibrator 3C147, which 
was also used as the bandpass calibrator. We spent ~95 minutes 
on target per epoch, or a total of ~ 668 minutes. 

The data were edited and calibrated using the Common 
Astronomy Software Applications (CASA 4.2.2) package, and 
the VLA calibration pipeline in its 1.3.1 version. After cal¬ 
ibration images were produced with pixel sizes of 0.2 arc- 
seconds, a natural weighting, and a multifrequency synthesis 
scheme (e.g.. Ran & CornwelH2011| l. The noise levels reached 
for each individual observation was about ~ 3 fiiy (see Table[^. 
Additionally, we produced images from the concatenated UV 
data from the epoch 2 to 7 and of the seven observations to reach 
lower noise levels of 1.02 /iJy and 0.95 /iJy, respectively. These 
noise levels are in agreement with the theoretically expected val¬ 
ues. 

Finally, we also performed the data reduction of three 
archived S-band (2-4 GHz) data sets taken with the VLA in its 
A configuration; these are part of the project 12B-06L Each in¬ 
dividual epoch runs for 1.0 hours. The receiver was also used 
in semicontinuum mode with the 8 bit sampler and 16 differ¬ 
ent spectral windows (with bandwidths of 125 MHz each) were 
recorded simultaneously to cover the full band. These observa¬ 
tions used the quasar 3C48 as the flux and bandpass calibrator, 
and quasar 12202-1-4216 as the phase calibrator. 

The data were edited, calibrated, and imaged following the 
same scheme as the new observations. The resulting sensitivity 
limits are ~ 11 //Jy and 6.6 /tJy, for the maps of individual epochs 
and for the map of the concatenated epochs (see also Table [^, 


' https://science.nrao.edu/facilities/vla/docs/manuals/obsguide 
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Fig. 1. VLA images with radio emission at the MWC 656 position. Left: image of the concatenated X-band data from Epoch 1 to 
7. Center: Epoch 1. Right: image of the concatenated X-band data from Epoch 2 to 6. Contour levels are -3, 2.8, 3.2, 3.8, and 4.2 
times the noise levels of the image (see Table[T]|. The red ellipse at the bottom left is the corresponding synthesized beam. The blue 
cross indicates the optical position of MWC 656. 


respectively. These values are also in agreement with those the¬ 
oretically expected. 


3. Results 


A source with a peak 3.4 times the noise level was detected in 
the first observation at orbital phase O = 0.49, i.e., at apastron 
passage. The source is coincident with the position of MWC 656 
(see Eig. center). A Gaussian fit to the source (using imfit in 
CASA) yields the flux density of the source to be 14.2 + 2.9 pJy. 
In the remaining observations, however, we do not detect any 
peak above 2.6 times the noise level of the images (see Eig. 
1^. By adding these six remaining observations, we detected a 
source with a peak 3.4 times the noise level, which is coinci¬ 
dent in position with MWC 656 and with the source detected in 
the first epoch (see Eig. [T] right). The flux density in this case 
was 3.7 +1.4 pjy. The two images of Eig. 1 (center, right) were 
produced with independent data sets, and support our detection 
of the radio counterpart of MWC 656. Interestingly, the flux in 
the first image (Eig. 1, center) has a flux density almost three 
times higher than that in Eig. 1, right, and this increased flux is 
occurring at apastron. We produce a final map by concatenat¬ 
ing the data of the seven observed epochs. The source is now 
detected at levels of 4.4 times the noise in the image (Eig. 1, 
left). The Gaussian fit to the source yields the flux density of the 
source to be 4.5 + 1.2 /liy, at the position RA=22^42™57.*305 
+0?005, DEC=+44°43'18'.'23+0!'08, and is consistent with a 
point like structure and parameters of Table 1. The source posi¬ 
tion is in good agreement with the optical position of MWC 656, 
RA=22''42™57.^30295, DEC=+44°43'1872525 (|van Leeuwen 
[2007l l. 


Einally, concerning the archived data (Table 1), we did not 
detect the source in any single S-band epoch, nor in the final 
image of concatenated observations at levels above 20 pJy. 




R.A. (2000.0) 



Epoch 7 

+ 


o 

_^^^ ^ 

R.A. (2000.0) 


Fig. 2. VLA images from epoch 2 to 7. The red circle at bottom- 
left is the synthesized primary beam. Contours are at -3, 2.6, and 
3 times the noise level. The blue cross at the center represents the 
position of MWC 656. The source remains undetected in each 
single image. The rms of the images and upper limits for the 
flux density of MWC 656 are given in Table 1. 
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4. MWC 656 in the context of the Lx, Lr relationship 


The radio flux values of 3.7+1.4 /rJy and 14.2 + 2.9 juJy at 
10 GHz, for a distance of 2.6 kpc, and assuming a nearly flat 
spectrum, corresponds to a Lr a t 8.6 GHz of 2.6x10^^ erg s~^ and 


( |2014| l de- 
‘ from an 


9.9 X 10^® erg s“\ respectively. |Munar-Adrover et al. 
termined an X-ray luminosity of Lx = 1-2 x 10'^* erg s 
observation at <1) = 0.08. We report our new data of MWC 656 
along with the measurements from Corbel et al. ( 2013| ). Figure 3 
shows (square/black) the radio and X-ray luminosities for the 
24 Galactic accreting binary BHs in the hard and quiescence 
states reported in Fig. 9 of Corbel et al. ( 2013| l. The position 
of MWC 656, given in red color, is rather close to the position of 
A062000 and XTE J1118+480, which are the weakest quiescent 
black holes known so far ( [Gallo et al.|2006||2014| l. 
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Fig. 3. Radio (8.6 GHz) vs X-ray luminosity (1-10 keV) dia¬ 
gram. Black squares indicate the positions of the 24 Galactic 
accreting binary black holes in the hard and quiescence states, 
as in Fig. 9 of Corbel et al. (2013; i.e., without the upper lim¬ 
its and neutron stars there present). The position of MWC 656 
is indicated by two red squares for the two fluxes of 3.7 //Jy and 
14.2 juJy. The position of MWC 656 is rather close to the faintest 
accreting black hole systems known so far, A0620-00 (triangle, 
Gallo et al.||2006|) and XTE Ill 18+480 (asterisk, [Gallo et al.j 
M4^ 


to test for MWC 656 the X-ray/radio correlation for black holes. 
The position of MWC 656 is close to the faintest black holes 
known so far, which are A0620-00 and XTE J1118+480. The 
secondary star in A0620-00 is a 0.7 Mq K3-K4V star in a 7.75 h 
orbit around the black hole ( [Gallo et al.|2006| l. The companion 
star in XTE J1118+480 is a K5-K8V star with an orbital period 
of 4.08 h ( [Gallo et al.|2014[ and references therein). MWC 656 
has an orbit of ~ 60 days and the companion star is a 10-16 
Mq Be star. Our observations support, therefore, the universal¬ 
ity of the Lx, Lr relationship, which is intimately related to the 
accretion-ejection coupling process, which seems to be invariant 
to different forms of accretion. 
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5. Conclusions 

Deep VLA observations of the BH-Be system MWC 656 have 
provided us with the first radio detection of an accreting stel¬ 
lar mass black hole with a Be star as companion. Along with 
a level of emission of 3.7+1.4 //Jy, determined combining six 
observations at the different orbital phases <1) =0.58-0.72 and 
fl) = 0.15-0.30, we detected in a single observation around apas- 
tron (<1) = 0.49) the flux density of 14.2 + 2.9 juJy. 


Several works have placed significant constrains on the 
Lx-Lr luminosities of quiescent systems (e.g., [Gallo et ah] 


[2006[ [20141 [Calvelo et al.[[2010| [Miller-Jones et al.|201 

1 . The 

nonsimultaneous XMM-Newton X-ray observation by 

Munar- 


[Adrover et al.[(|20T4|) at orbital phase O = 0.08, also allowed us 
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